General Cramer-Rao bound for parameter estimation using Gaussian multimode 

quantum resources 
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Multimode Gaussian quantum light, including multimode squeezed and/or multipartite quadra- 
ture entangled light, is a very general and powerful quantum resource with promising applications 
to quantum information processing and metrology involving continuous variables. In this paper, 
we determine the ultimate sensitivity in the estimation of any parameter when the information 
about this parameter is encoded in such Gaussian light, irrespective of the exact information ex- 
traction protocol used in the estimation. We then show that, for a given set of available quantum 
resources, the most economical way to maximize the sensitivity is to put the most squeezed state 
available in a well-defined light mode. This implies that it is not possible to take advantage of the 
existence of squeezed fluctuations in other modes, nor of quantum correlations and entanglement 
between different modes. We show that an appropriate homodyne detection scheme allows us to 
reach this Cramer-Rao bound. We apply finally these considerations to the problem of optimal 
phase estimation using interferometric techniques. 

PACS numbers: 03.65.Ta, 42. 50. Ex, 42.50.Lc, 42.50. St 



Optical techniques are widely used in many areas of 
science and technology to perform high performance mea- 
surements and diagnostics, for example in spectroscopy, 
ranging, trace detection, microscopy, image processing. 
There are many reasons for this: light allows us to extract 
information in a remote and non destructive way, it car- 
ries information in a massively parallel way, and perhaps 
more importantly optical measurements can reach very 
high precision and/or sensitivity levels, the best exam- 
ple being the interferometers used as gravitational wave 
antennas. In this respect, it is very important to know 
exactly and in the most general way what are the limits of 
accuracy in parameter estimation using light. These lim- 
its depend on the noise present in the detection process, 
which is ultimately, when all sources of imperfections in 
the setup have been eliminated, the quantum fluctua- 
tions of the photodetection signals used to estimate the 
parameter of interest. 

When the detected light is produced by ordinary, shot 
noise limited, sources, such a limit is called "standard 
quantum limit" . It is well-known that it is possible 
to improve optical measurements beyond the standard 
quantum limit by using squeezed [l[ or entangled light. 
This statement has been demonstrated for what can be 
called " simple measurements" , in which the information 
about theparameter p to measure is carried by the total 
intensity^ or by the phase[H, [B| of a single-mode light 
beam. It has been recently extended to the case of opti- 
cal images in which the parameter p to measure does not 
change the total intensity of the light but instead modi- 
fies the details of the repartition of light in the transverse 
plane0, and experiments have been performed for small 
displacement measurements Q. 



The purpose of this paper is to derive the expression of 
the ultimate sensitivity in the estimation of a parameter 
encoded in any kind of variation, spatial and/or tempo- 
ral, of an optical field. To this purpose, we will deter- 
mine the Cramer-Rao bound on the parameter estima- 
tion, which gives a limit to the best expectable precision, 
independently of the strategy used in the estimation 0- 
Hoj . In addition, we will consider that the light which is 
used in the measurement belongs to a wide set of quan- 
tum states, namely the set of multimode Gaussian states, 
which contains all possibilities of multimode squeezing 
as well as all kinds of multipartite quadrature entangle- 
ment. These non-classical states are now produced by 
experimentalists with impressive amounts of squeezing 
or entanglement [ll[. We will also assume that the in- 
formation that is extracted from the light and then pro- 
cessed is the intensity, or a field quadrature, at different 
positions and/or times. Our approach does not giv e the 
most general " quantum Cramer- Rao bound" [12l4l4j , op- 
timized over all possible quantum states and all possible 
quantum measurements, but instead the lowest possible 
uncertainty optimized over a subset of quantum states 
and quantum measurements that are readily available us- 
ing present technology even for very large values of the 
mean photon number (up to 10 16 ), which is not the case 
with more exotic non-Gaussian states and detectors like 
NOON states and Photon Number Resolving detectors. 

A sketch of a general optical measurement setup is dis- 
played in the figurcQ] An optical device, containing light 
sources and various optical elements, produces an output 
beam of light described at the detector position by its 
complex electric field E ou t, which is multimode in space 
and time in the general case. The parameter that one 
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FIG. 1: General scheme for estimating light parameters 



wants to determine, called p, modifies in a known way 
the output beam, so that E out = E out (p). One then 
uses a photo-detector that records a p-dependent infor- 
mation. This information is a classical quantity which is 
processed by analog or digital techniques. The outcome 
is a real number which is used as an estimator of p, that 
we assume to be unbiased. 

In the following, we assume that both the optical de- 
vice and the photodetector are perfect, i.e. that all causes 
of technical noise have been eliminated. The only source 
of uncertainty in such a measurement is thus the un- 
avoidable fluctuations in the photodetection signal aris- 
ing from the quantum fluctuations of the light beam E out , 
and that can be tailored in the optical device producing 
the output beam in such a way that E out may have either 
squeezed or quantum-correlated quadrature fluctuations, 
or both. 

To simplify the formalism we will assume that the out- 
put field is polarized in a fixed direction, so that one 
can use a scalar description of the output field E out (p). 
In addition we will use the paraxial and slowly varying 
envelope approximations when necessary. We will also 
assume that the number of modes necessary to describe 
the output field is finite, and has a value M. The mean 
output complex field can then be written 



E ont (r,t,p) = i 



Nfiw 
2e cT 



ui(r,t,p) 



(1) 



where u\(r^t,p) is a normalized mode that follows the 
evolution of the mean field with p, T the exposure time 
and N the mean value of the total photon number. We 
will leave aside the already well known case where the 
parameter is encoded in the light intensity, and assume 
that N is not modified by the parameter p. 

We introduce a general mode basis to describe the 
detection process that we call {vi(r, t)} (i = 1, ...,M), 
which constitutes a complete basis for the description 
of the electric field at the location of the photodetector. 
One can write the positive frequency electric field oper- 
ator 



hiij 
2e a cT 



(2) 



where hi is the annihilation of photons operator in mode 
V{. The modes Vi(r,t) can be for example the Hermitc- 



Gauss modes in the pure spatial domain, or a set of Gaus- 
sian temporal modes if one wants to describe light pulses 
of arbitrary shapes. We write \^{p)) the quantum state 
produced by the device and which depends on p. One has 
E out {r,t,p) = (ilj(p)\E^{{r,t)\tl){p)) (the bar denoting in 
all this paper the quantum mean value). We also intro- 
duce the quadrature operators for mode i: Y L + = h\ + hi 
and Yf = i(b\ — hi) and the general matrix column 
quadrature operator Y = (Y^~ , Y^, Y{~ , Y^) T . 

We assume that \ijj(p)) is a Gaussian state of light. It 
can then be described by a Gaussian Wigncr function 
and its properties are entirely defined by the mean value 
and covariance matrix of its quadratures operators. Let 
us call Y(p) the mean value of the quadratures and S(p) 
the covariance matrix written in the basis {v{, 6^}, which 
contains all the information about the multimode entan- 
glement and squeezing properties of the state. Its Wigner 
function is then expressed as 



Wp(Y) 



(27T) 



M 



,--(Y-Y(p)) T S- 1 (p)(Y-Y(p)) 

(3) 

In order to get simple expressions, we assume that we 
have a pure state (det[S] = 1), but one can easily gener- 
alize the derivation of the Cramer-Rao bound to a mixed 
state. 

Let us now calculate the Fisher information, inverse of 
the Cramer-Rao bound, corresponding to the detection 
of a variation of the parameter p from its initial value, 
that one can always take as p = by a change of ori- 
gin. We need to assume here for the calculation that one 
is able to measure the two quadratures of all the modes 
of the field. Even if non-physical, as two quadratures 
of the same mode correspond to non commuting observ- 
ables, this assumption will give an upper bound to the 
Fisher information as it increases the information that 
one can actually extract from the field. More physically, 
for a stationary situation where the same state of light 
is continuously produced by the device, it corresponds 
to measurements performed during twice the time, once 
on each quadrature. This means that we have complete 
access to the vector Y. In this case, and assuming that 
we have only unbiased estimators^ the Fisher information 
has the following expression jT 



-^Fisher 



d 2 l(Y\p) 



dp 1 



L(Y\p = 0) alY 



L{Y\p) being the likelihood of measuring Y knowing 
that the parameter value is p and l{Y\p) being the log- 
likelihood. 

Using the Wigncr function defined in eq. [3] one finds 
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(4) 



This expression gives the Fisher information as a function 
of the field mean value and variance, and of its depen- 
dence with p. It is made of two terms. The first term is 
linked to the variation of the mean field distribution with 
p while the noise taken into account is the one for p = 0. 
The second term comes from the modification of the co- 
variance matrix of the field with p in the neighborhood 
of 0, and it does not depend on the mean field, thus it 
does not scale with the mean number of photons N. For 
large value of N, as those involved in continuous variable 
regime, it becomes negligible compared to the first term 
of the Fisher information, and we will neglect it in the 
rest of this paper. In the remainder of this paper, we will 
also restrict our analysis to the case of the estimation of 
very small variations of the parameters value, where the 
effort towards precision estimation enhancement is most 
crucial. 

We now choose an adequate detection basis in order to 
obtain an interpretable form of the Cramer-Rao bound. 
The idea is to find a basis where ^p- 1 „ is single mode. 

dp lp=0 ° 

Let us define the first mode of the detection basis by 



Vi = p c 



dui(p) 



dp 



where — 

P =o P c 



dui(p) 




dp 


p=0 



(5) 



and put no additional constraints on the other modes 
v n >\ of the basis (except that they are orthogonal to 
v\). Let us stress that this basis does not depend on the 
actual value of p. One has ^| p _ = 2V~Nvi/p c . The 
calculation of the Fisher information in that basis is then 
very simple as only one vector is involved and thus only 
one element of S _1 (0). 



(6) 



where E^^O) is the first left, top element of the matrix 
Y,~ l (p = 0). We find that the Fisher information only 
depends on the mode v\ defined by the derivative of the 
mean field mode with respect to p, that we will call from 
now on the detection mode [THJ] . 

We are then led to our final expression of the Cramer- 
Rao bound for unbiased parameter estimation using 
quantum Gaussian states: 
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1 



Sa!i)(0) 



(7) 



The Cramer-Rao bound we have just derived depends 
only on three parameters: a scaling factor p c that charac- 
terizes the variation of the mean field distribution with p, 
the mean total number of photons and one element of the 
inverse of the covariance matrix. In this article we take 
for given p c and N as optimizing the mean field shape 
and increasing the number of photons are two well known 
strategies to increase information extraction and we as- 
sume here that they have both been optimized. We focus 
here on the effect of the noise, and the remarkable result 
that we have obtained is that the Cramer-Rao bound de- 
pends only on the diagonal matrix element of the inverse 
of covariance matrix on the detection mode. 

The Cramer-Rao bound gives the best estimation 
precision one can hope to obtain on a parameter en- 
coded on a specific Gaussian quantum state, character- 
ized by its covariance matrix. The question is now to 
find the best strategy to optimize this bound given a 
certain amount of quantum resources. Let us start, like 
in most experiments to date, from s squeezed vacuum 
states, erj (i = 1, ...s) is the r.m.s. value of the squeezed 
quadrature of mode i, and we call er m i n the smallest of 
its values. With the help of linear couplers i.e. of uni- 
tary transformations [l6[ on the mode basis, the multi- 
mode squeezing can be transformed partially or totally 
into multipartite entanglement in a different mode basis. 
One can show that, under such unitary transformations, 
the diagonal coefficients of the inverse of the covariance 
matrix are bound by its spectral radius, which is equal to 
i/ CT min- The equality is reached only when the detection 
mode 1 is an eigenmode of the covariance matrix with 
the eigenvalue cr 2 ^, and thus when the detection mode is 
the most squeezed mode and is not correlated with all the 
other modes. The optimized Cramer-Rao bound is thus 



Ap 



P t 



Pc 



2VN 



(8) 



We have shown here that the only way to minimize 
the Cramer-Rao bound given a certain initial amount 
of quantum resources is to put the most squeezed state 
available in the detection mode and not to have correla- 
tions with other modes. The presence of other squeezed 
modes, or of any kind of entanglement, will not help 
to improve the limit: one cannot take advantage of 
squeezed fluctuations or quantum correlations coming 
from different modes to improve the estimation of a single 
parameter 17 1. 

The determination of the Cramer-Rao bound does not 
tell us how to reach it, or even whether it is actually 
possible to reach it. In the present case, our analysis has 
outlined the importance of what we have called the detec- 
tion mode. We will now show that a balanced homodyne 
detection scheme in which the local oscillator is taken in 
the detection mode allows us to reach the Cramer-Rao 
bound. 

When p is close enough to its initial value 0, the mean 
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field mode ui(j>) can be expanded as 



P 



= ui(0) + ^t>i (9) 

p=0 Pc 



For a small non-zero value of p, the field is therefore 
composed of the original field plus the detection mode 
Vi times the parameter. In the general case, v± is not 
orthogonal to Ui(0). We define a detection orthonormal 
basis whose first mode is v\ and second mode, in the sub- 
space spanned by {ui(0), v{\, is v 2 = ^gj^gfefep ■ 
Noting that J ul(0)^\ p —0 is an imaginary number, we 
can write ui(0) = icuVi + c\ 2 v 2 where en is a real num- 
ber. Finally, u\ (p) can be expressed in the detection basis 
as 



Mp) = ( — 



+ lCn V\ + C\ 2 V 2 



(10) 



On can show that[18|, using a homodyne detection 
scheme in which the local oscillator is a coherent state 
in the v\ mode, and if the relative phase between the 
field to detect and the local oscillator is zero, the inten- 
sity difference operator between the two outputs of the 
set-up is given by 



I- = 



Hid 

2e cT ' 



'N 



Pc 



(11) 



where Nq is the local oscillator mean number of photon 
and 6Y1 = — (Yi) is the quantum noise fluctuation 
operator of the incoming field in v\ mode. From that 
expression, it is easy to derive both the signal (7_) and 
the noise Uj . The p value giving a signal to noise ratio 
equal to 1 is 



PSQL 



2VN 



dui 



\p=0 



(12) 



which is nothing else than the Cramer-Rao bound (|8]| 

We will finally illustrate the interest of our approach by 
revisiting a well-known and important problem of quan- 
tum optics, namely the interferometric measurement of 
a phase shift <f> [19h22| . As a Michelson interferometer 



has two output modes v\ and v 2 than can be detected 
separately, it is a two-mode problem. In presence of a 
phase-shift <p between the two arms of the interferometer 
close to a bias value <f>o, the total mean output field is 



Ul {<j>) = Vi cos(F(</>/2)) + v 2 sin(F(</>/2)) 



(13) 



When the arms are empty, F(x) = x, but optical devices, 
such as Fabry-Perot cavities, can be inserted in the two 
paths in order to increase the phase sensitivity of the in- 
terferometer. A complete basis of our two-mode space is 
made of the mode u\ completed by the orthogonal mode 



u 2 (4>) = -uisin(F(0/2)) + v 2 cos(F((/>/2)). The detec- 
tion mode defined in ([5]) is in the present case u 2 ((f>o), 
and the Cramer- Rao bound, according to equation©, is 



1 



\F'\VWs- 



(14) 



It is independent of the initial phase bias 4>o , and is min- 
imum when F' is maximum, as could be expected. Its 
minimum value is obtained by having at the output of the 
interferometer a quantum state consisting of a squeezed 
state in the detection mode u 2 (4>o) and any Gaussian 
state in the orthogonal mode ui(<f>o) uncorrelated with 
the first one. It is easy to see that such an output state 
is obtained by sending at the two input ports of the inter- 
ferometer a tensor product of a field with nonzero mean 
value and of a vacuum squeezed state: one thus finds 
that the well-known scheme introduced by C. Caves Q is 
optimal when one uses Gaussian resources. This implies 
in particular that the use of Gaussian entangled states of 
the two input modes will not help improve the sensitivity 
of the interferometer 17 1 . 

Let us finally mention that these results can be gener- 
alized to the estimation of several parameters: one can 
show that when these parameters are " orthogonal" , it is 
possible to optimally reduce the noise on their estima- 
tors at once by independently squeezing the two corre- 
sponding detection modes. We now plan to extend the 
present approach to various kinds of non-Gaussian states 
and detectors in the continuous variable regime, such as 
Schrodinger cats and multiphoton detectors. 
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